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ABSTRACT 

Using Sloan Digital Sky Survey Data Release 5, we have investigated the intrinsic axis ratio dis- 
tribution (ARD) for early-type galaxies. We have constructed a volume-limited sample of 3,922 
visually-inspected early-type galaxies at 0.05 < z < 0.06 carefully considering sampling biases caused 
by the galaxy isophotal size and luminosity. We attempt to de-project the observed ARD into three- 
dimensional types (oblate, prolate, and triaxial), which are classified in terms of triaxiality. We 
confirm that no linear combination of randoms-distributed axis ratios of the three types can repro- 
duce the observed ARD. However, using Gaussian intrinsic distributions, we have found reasonable 
fits to the data with preferred mean axis ratios for oblate, prolate, and triaxial (triaxials in two axis 
ratios), pt = 0.44, pt v = 0.72, ^ t ,f3 = 0.92, fi t ,-y = 0.78 where the fractions of oblate, prolate and triaxial 
types are 0:P:T = 0.29 ±0 09 : 0.26 ±o n : 6.45 ±013 . We have also found that the luminous sample 
(—23.3 < M r < —21.2) tends to have more triaxials than the less luminous (—21.2 < M r < —19.3) 
sample does. Oblate is relatively more abundant among the less luminous galaxies. Interestingly, the 
preferences of axis ratios for triaxial types in the two luminosity classes are remarkably similar. We 
have not found any significant influence of the local galaxy number density on ARD. We show that 
the results can be seriously affected by the details in the data selection and type classification scheme. 
Caveats and implications on galaxy formation are discussed. 

Subject headings: surveys — method: statistical — galaxies: elliptical and lenticular, cD — galaxies: 
formation — galaxies: kinematics and dynamics 



1. INTRODUCTION 

Ever since Hubble (1926) investigated the apparent 
flattenings of early-type galaxies, numerous studies have 
attempted to pin down their intrinsic shapes. Such ef- 
forts mainly focused on the apparent axis ratio not only 
because it is easy to measure but also its distribution can 
be used to extract the kinematics and may even constrain 
the general formation history of galaxies. 

Early studies using apparent axis ratio distribution 
(ARD) were made on the assumption that early-type 
galaxies are composed of one type. Assuming oblate- 
ness alone, Sandage, Freeman, & Stokes (1970) suggested 
the possibility of Gaussian or skewed binomial distribu- 
tion of intrinsic ARD using Reference Catalog of Bright 
Galaxies (RC1, de Vaucouleurs & de Vaucouleurs 1964). 
Binney (1978) 's approach assuming prolateness also suc- 
cessfully reproduced the apparent ARD from RC1. The 
existence of triaxial galaxies (Bertola & Capaccioli 1975; 
Illingworth 1977) opened up a better chance to repro- 
duce the observed ARD. Binggeli (1980) and Benacchio 
& Bertola (1980) showed that the apparent ARDs for 160 
ellipticals from Revised-Shapley-Ames catalog of bright 
galaxies (Sandage & Tammann 1979) and 348 ellipticals 
from Strom & Strom ( 1978a, b,c) were well represented 
as a group of triaxial galaxies assuming a fixed case of 
(3 = (1 + 7)/2 where (3 is the ratio between the second 
longest to the longest axis and 7 is the shortest to the 
longest axis. In particular, Binney & de Vaucouleurs 
(1981) aimed to reconstruct the apparent ARD for RC2 
(de Vaucouleurs et al. 1976) using Lucy's inversion tech- 
nique (Lucy 1974, see also Noerdinger 1979). They con- 
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sidered oblate, prolate and triaxial (0,P,T) separately 
and still found acceptable fits to the observed ARD. The 
limitation of this technique, however, is that it is difficult 
to constrain the ratios between O, P, T galaxies. 

More recently, Fasano & Vio (1991) concluded that a 
purely-biaxial model cannot reproduce the small num- 
ber of apparently-round galaxies. This paucity however 
looks significantly different when samples are drawn from 
different catalogs (e.g., RC1, RC2, or Revised-Shapley- 
Ames Catalog of Galaxies). This obviously results in 
disparate intrinsic distributions. Furthermore, Lambas 
et al. (1992) found a reasonable fit to the sample of 2,135 
galaxies from the APM Bright Galaxy Survey (Maddox 
1990) assuming all ellipticals are triaxial whose intrinsic 
distributions are two dimensional Gaussian. The similar 
work done by Ryden (1992) presents consistent results. 
These results support the assertion that the paucity 
of round galaxies can be reproduced by a dominantly- 
triaxial galaxy population. Considering its usefulness, it 
is important to accurately sample the observed ARD. 

Fortunately, Sloan Digital Sky Survey Data (Adelman- 
McCarthy et al. 2005) allows us to study the apparent 
flattenings for a large number of galaxies. Our effort 
to make a complete volume-limited sample is one of the 
main factors distinguishing our work from previous stud- 
ies. 

We also note the interesting result of Davies et al. 
(1983) that bright ellipticals may be more slowly rotat- 
ing than faint ones. Tremblay & Merritt (1996) found 
that the apparent ARD is markedly different for two lu- 
minosity classes ; the ellipticals brighter than Mb — — 20 
are rounder than less luminous ones. In particular, re- 
cent observations show that the division clearly occurs at 
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Mb = —20.5 (Rest & van den Bosch 2001) and also sup- 
port previous works suggesting that bright galaxies show 
a core profile and boxy isophotes while faint galaxies have 
a power-law profile and disky isophotes (Bender 1988, 
Kormendy & Bender 1996, Faber et al 1997). Therefore 
we investigate on the effect of the galaxy luminosity on 
the ARD. 

Another factor we should note for intrinsic shape 
of early-type galaxies is an environmental dependence. 
Dressier (1980) pointed out a density-morphology rela- 
tion which reflects the importance on the formation pro- 
cess. We thus search for the connection between intrinsic 
shape and environment for early types. 

In this paper, we propose two simplifying assumptions; 

(i) early-type galaxies are geometrically perfect ellipsoid, 

(ii) they are randomly oriented. We assume that early 
types consist of oblate, prolate and triaxial. We believe 
this assumption makes our approach more realistic than 
the previous models composed of only one or limited 
types. On this basis we investigate the projection effect 
on the apparent ARD. We simply describe the model 
distribution for oblate, prolate and triaxial in §2. We in- 
troduce our sample selection with completeness tests in 
§3. In §4, we investigate the intrinsic ARD for volume- 
limited samples. In §5, we analyze the intrinsic shape of 
two different luminosity samples. Dependence of envi- 
ronment on ARD is investigated in §6. In §7, we discuss 
the limitations of our approach. Finally, we discuss re- 
sults and their implication in §8. 

2. ANALYTIC APPARENT AXIS RATIO DISTRIBUTION 

We use an analytical description to calculate the prob- 
ability distribution (Franx et al. 1991; Binney & Mer- 
rifield, 1998). The probability of finding the apparent 
ellipticity (e) in the interval (e, e + de) is 

' J' 11 dfxde for < e < e u 

for ei < e < 1, 

(1) 

where e = (1 - e) 2 , h(r) = (r - a 2 )(r - b 2 )(r - c 2 ) 
and ei = 1 — b/a. a, b and c indicate three axes of 
an ellipsoid, respectively. For galaxies with b/a > c/b 
(oblate-triaxial), /ii and /i 2 , which depend on the galaxy 
shape, are 

[ (b 2 ,b 2 /e) for 0<e<e 2 , 
(Mi,M2) = < {b 2 ,a 2 ) for e 2 < e < e 3 , (2) 
[ (c 2 /e, a 2 ) for e 3 < e < t\ 

where £3 = 1 — c/b. For galaxies with b/a < c/b (prolate- 
triaxial), /xi and /i 2 can be written as 

( (b 2 ,b 2 /e) for 0< e < e 3 , 
Oi., ^2) = < (c 2 /e,6 2 /e) for e 3 < e < e 2 , (3) 
[ (c 2 /e,a 2 ) for e 2 < e < a 

In order to illustrate the type (OPT)-dependence of the 
apparent ARD, we adopt the classification scheme of 
Franx et al. (1991). To classify the early-type systems, 
we use triaxiality (T), 

1 - (3 2 
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and each type can be expressed as 
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Fig. 1. — Scheme of galaxy classification. In this study, we adopt 
Franx(1991)'s scheme. Shaded area : Based on the observational 
constraints, we use axis ratios greater than 0.2. 




Fig. 2. — Apparent axis ratio distribution, <E> (i.e., the probability 
density of the projected shapes of a group of galaxies) for each type 
assuming a uniform intrinsic ARD. Solid, dotted, and dashed lines 
correspond to oblate, prolate, and triaxial, respectively. 



In this study, we assume that there is no early-type 
galaxy with axis ratios smaller than 0.2 because such 
systems are rarely observed. In Fig. [T] we display the 
classification scheme. Fig. [2] shows the simplest special 
case of the apparent ARD, that is based on (uniformly- 
distributed intrinsic axis ratios. We mean no preferred 
values of the intrinsic ratios by "uniformly- distributed" . 
The numbers of the samples of three (OPT) types simu- 
lated are normalised to be the same, hence unbiased by 
the area difference between the types in FigfT] If an ob- 
served ARD has a large number of round galaxies near 
1, we can deduce that oblate is the main component, for 
example. 

3. SDSS SAMPLE SELECTION AND DATA ANALYSIS 

The SDSS provides a large homogeneous database. 
There are approximately 28,000 galaxies within 0.05 < 
z < 0.06, an excellent sample for studying ARD statis- 
tically. In this section, we describe our data selection 
scheme and completeness tests. 

3.1. Morphological Classification 
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A well-defined criterion for morphological classification 
is necessary to study the ARD of early- type galaxies. Be- 
cause visual inspection of all galaxies is extremely time- 
consuming and still subjective, we adopted the SDSS 
pipeline parameter fracDev, which indicates the frac- 
tion of the brightness profile that can be explained by 
the de Vaucouleurs profile (de Vaucouleurs 1948) as fol- 
lows: 

f composite = fracDevf deV + (1 - fracDev)f exp (5) 

where fd e v indicates de Vaucouleurs fluxes. We assume 
that a conservative sample of early-type galaxies have 
fracDev > 0.95 in all 3 bands, g'r'i', following the prac- 
tice of Yi et al. (2005). Using this fracDev parameter, 
we compile 4,994 galaxies within 0.05 < z < 0.06. 

3.2. Data analysis 
3.2.1. Luminosity Dependence 

A complete volume-limited sample is crucial for this 
study because it has a direct effect on the intrinsic shapes 
of galaxies. But SDSS provides spectroscopic informa- 
tion only for the galaxies of r < 17.77 ; hence, our sam- 
ple cannot be free from luminosity bias. To investigate 
this, we need to know how the ARD varies with the size 
of major axis. Fig. [3] shows the trend that the luminos- 
ity gradually increases with increasing minor axis for a 
fixed major axis size. This effect is clearer for the larger 
galaxies (red and blue dots) but less clear for smaller 
galaxies (black dots) probably because of the magnitude 
limit. We also denote the general trends of the three 
different luminosity classes with three lines. It is clear 
that small faint flat early types are more easily missed, 
and so the luminosity limit biases the apparent ARD. 
We have decided to exclude galaxies with major axis ra- 
dius IsoA r smaller than 16.2". This however has a ten- 
dency of removing distant faint galaxies from the sam- 
ple. In order to alleviate this problem, we construct a 
volume-limited sample by selecting close galaxies within 
0.05 < z < 0.06. Obviously, no redshift dependence 
is considered important here. Note that Odewahn et al. 
(1997) already reported that there is no significant differ- 
ence in the apparent ARD between distant and nearby 
samples. Constraining the absolute magnitude range, 
statistical Kolmogorov-Smirnov (KS) test confirms that 
distant (0.097 < z < 0.1) and close (0.05 < z < 0.06) 
galaxies share the same parent ARD with a 99%-level 
confidence. Therefore, the results from the analysis on 
our close (0.05 < z < 0.06) sample likely holds for a 
larger redshift range. 

3.2.2. fracDev Dependence 

It is also of interest whether our main selection crite- 
rion fracDev > 0.95 influences the ARD. S. Joo (priv. 
comm.) performed an independent classification of early 
types using fracDev and visual inspection, for two dif- 
ferent fracDev limits, fracDev > 0.95 and 0.50 <fracDev 
< 0.95. He pointed out to us that the higher fracDev cri- 
terion misses some early types, mainly flatter galaxies. 
For example, if one uses lower Fracdev limit, the fraction 
of relatively flatter (b/a < 0.6) galaxies with M r > —20.5 
would be 5.9%, while the fraction for higher fracDev limit 
is 3.9%. In this regard, our sample does not represent the 
entire early-type galaxy population but is slightly biased 
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Fig. 3. — The sample bias caused by the optical brightness 
criterion. Gray dots in the background indicate all galaxies within 
0.05 < z < 0.06. Galaxies of different apparent isophotal major- 
axis sizes are grouped into 3 bins: orange dots denote the largest 
galaxies. Three dashed lines show the trend that rounder galaxies 
are brighter for a fixed apparent axis ratio. Faint flat small galaxies 
(IsoA r < 16.2"), near black dotted line, are missed in this kind of 
survey (see text). 

towards rounder galaxies. We decided, however, not to 
worry about this, first because the exact fracDev crite- 
rion for early types is unclear, and second because we are 
for the moment more interested in the methodology. 

3.3. Test Sample 

Our final sample of 3,922 galaxies is chosen with red- 
shift criterion (0.05 < z < 0.06) and major axis criterion 
(IsoA r > 16.16"). We remove a small number of rel- 
atively faint outliers of r > 17.5 to ensure reasonable 
image quality. But this has no impact on our results 
at all. We also exclude 1,072 galaxies from the 4,994 
galaxy sample because they appear to be spiral contami- 
nants or severely-distorted in the visual inspection. The 
ARD for final sample is shown in Fig. [4] For the pur- 
pose of comparison, we also plot the apparent ARD from 
APMBGS data (Loveday 1996). We bin the data by the 
size roughly drawn from the Izenman method (1991), 
alQRn^ 1 / 3 , based on the total number n and interquar- 
tile range (IQR) where 2.0 < a < 2.5 (Izenman 1991). 
The peak around p — 0.8 is noteworthy. 

4. INTRINSIC AXIS RATIO DISTRIBUTION FOR 
VOLUME-LIMITED SAMPLE 

It was demonstrated in Fig. [5] that the apparent ARD 
is not uniform even though the intrinsic distribution is 
uniform. To extract the intrinsic shapes of our galax- 
ies, we use composite models of O, P, T types based on 
the Franx et al. (1991) classification scheme (Fig. QJ. 
In this section, we compare the observation with model 
distributions using two strategies. 

4.1. Uniform distribution of intrinsic axis ratio 

We investigate the possibility of uniform intrinsic dis- 
tribution. To measure the goodness of the fit, we use re- 
duced x 2 , and this is constructed by considering Poisson 
error. Using the volume-limited sample from §3, we try 
to derive the fractions of the O, P, T types in the combi- 
nation that best reproduces the observed ARD. We can 
express this with weight (IF) of each type. 
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Fig. 4. — Our final SDSS sample compared to the previous APM 
sample of Loveday (1996) (Gray dots). We derive a relatively- 
complete sample concerning luminosity and major axis radius. 



<f> W + %W P + §tW t = <P obs (6) 

where apparent ARD = 7) 

7 

Index i indicates each type, and \& is the apparent ARD 
for certain axis ratio (,3,7). Weight directly reflects the 
frequency of each type. The best solution that makes a 
minimum X 2 is W Q : W p : W t = 0.09 : 0.00 : 0.91 but 
for a very poor statistic (x 2 ed ~ 50). Hence, no linear 
combination of OPT types with random ARD reproduces 
the observed ARD of our sample!. 

4.2. Gaussian distribution of intrinsic axis ratio 

Since the random intrinsic ARD fails to reproduce 
the observed data, we adopt Gaussian distribution. Al- 
though previous studies, which allow arbitrary distribu- 
tion, could estimate the preference of axis ratio, it is 
nearly impossible to quantify the fraction of oblate, pro- 
late and triaxial type. In this respect, Gaussian is the 
best distribution to test the fraction and perference of 
axis ratios. The Gaussian distribution can be written as 

1 (x — [l) 2 
Fgau(x; fi, a) = —^=- exp [ ^T~\ ( 7 ) 

where fj,, a correspond to the mean position and width of 
Gaussian distribution. Then preferences for oblate and 
prolate galaxies can be expressed as 

$o = Y.Y. F 9<™(r,Vo,cr )<f ((3,-f) (8) 

7 

7 

For triaixal galaxies, we use two Gaussian weights fol- 
lowing the Lambas et al. (1992) approach. 

7 

(9) 

for {08,7) I T ijmin < (±5^) < T hmax } where 0.2 < 7 < 
< 1.0. 



To reproduce the observed ARD above, we imposed 
specific intrinsic axis ratios and sum over three probabil- 
ity distributions after multiplying their weight factors. 
We should note that these weight factors are different 
from Gaussian weight in Eqn. [6j and recall that this 
provides a simple framework where we can investigate 
the fraction of each type. 

With this approach we find good matches. The 
best-fit model that yields minimum x 2 ~ 1 is 
(jj. , lip, li t ,0, Mt )7 ) = (0-46, 0.72, 0.92, 0.74), and 
(a a , a p , a t ,p, 04,7) = (0.1, 0.05, 0.1, 0.3). For this 
case, we found that the fraction of each type is 0:P:T = 
0.35 : 0.18 : 0.47. However, our parameter space is so 
complicated that the minimum \ 2 model may not repre- 
sent the most meaningful result. Instead, the statistical 
properties of all possible models with reduced \ 2 values 
within Id range ( A% 2 = xt ~Xmin — -0 are more mean- 
ingful because they all show good agreements with the 
observation (Fig. [5]). Our results show that the triaxial 
component is dominant around the high axis ratio, while 
oblate also plays an important role in the low axis ratio 
region. Statistically, the total fraction of each type is 
0:P:T = 0.29 ±009 : 0.26 ±o n : 0.45 ±013 in 1 a range, 
which suggests that triaxial early-types are most com- 
mon. For comparison to the "best-fit model", the % 2 
space of the "good models" (A% 2 < 1) show a conver- 
gence at a slightly different configuration (/i D , lit 0, 
= (0.44, 0.72, 0.92, 0.78) as shown in Fig. [5] We 
believe that this is a more statistically-representative re- 
sult. 

In the similar simulation of Lambas et al. (1992), 
the optimal solution for fitting the apparent ARD 
from APMBGS data with two-dimensional Gaussian had 
l^t,p = 0.95, nt.-y = 0.55, <7t,0 = 0.35, <7 tn — 0.2. Likely 
reasons for the difference between their results and ours 
are: (1) they considered only triaxial and used differ- 
ent classification scheme; and (2) the observed ARDs are 
slightly different. They used smaller values of /Zt i7 (i.e., 
flatter) than ours probably in order to fit the low axis 
ratio regions without considering oblate and prolate ele- 
ments. On the other hand, Ryden's (1992) results that 
produce the best-fit model with fi t ^ = 0.98, /x t:7 = 0.69 
and at = 0.11 are closer to our results; but, since their 
observed ARD was suppressed in the low axis ratio re- 
gion, their best-fit triaxials were rounder than those in 
the Lambas et al. fit. 



5. INTRINSIC ARD FOR DIFFERENT LUMINOSITIES 

We investigate whether two luminosity classes men- 
tioned in §1 (Bender 1988, Kormendy & Bender 1996) 
have different intrinsic shapes by comparing their best-fit 
(within 1 a) model OPT ratios. First, we divide galax- 
ies into "luminous" (M r < —21.2) and "less luminous" 
{M r > -21.2)) groups. Rest & van den Bosch (2001) 
found a division at Mg ~ —20, which corresponds to 
M r ~ —21.2 using the transformation of Smith et al. 
(2002) and the typical color for early-types (B—V ~ 0.9). 
As seen in Fig. [71 the "luminous" sample exhibits a 
larger number of round galaxies. We derive the weight 
and preference of intrinsic axis ratio for each type. 

We first focus on the weight of each type between the 
two groups. Averaged over the region of the la range 
(A% 2 < 1), the OPT weights for the "luminous" and 
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Fig. 5. — Apparent axis ratio distributions based on our Gaussian 
weight scheme. The OPT-combined best fit is shown by the dashed 
line, and the corresponding component fits (O, P, T) are shown. 
We exhibit the average value with the dotted line and 1<t deviation 
of each model that satisfies Ax 2 < 1 as shaded region. Note that 
triaxial galaxies are the dominant component. 

"less luminous" samples are as follows. 

O:P:T = 0.13 ±008 : 0.20 ±013 : 0.67 ±013 (luminous) 

O:P:T = 0.38 ±008 : 0.18 ±01 ° : 0.43 ±011 (less luminous) 

This implies that luminous early-types are likely triaxial, 
while there still exists a large amount of oblate galaxies 
in the "less luminous" sample. We display in Fig. [8] 
how each type can be viewed in the sky. Note that the 
galaxies with a high apparent axis ratio are likely tri- 
axial regardless of their luminosity, and oblate galaxies 
are more common in the "less luminous" sample. This 
may indicate that different formation process between 
the two samples. We will discuss this in greater detail 
in §8. The dichotomy, if real, might be explained by the 
argument presented by Valluri & Merritt (1998) involv- 
ing the central supermassive black hole and the crossing 
time difference between the bright and faint ellipticals. 

However, regarding preferences of intrinsic axis ratio 
(/z) and Gaussian widths (a), there is no clear disparity 
between the two samples. The "luminous" sample shows 
a preference for /it i( g = 0.90, /it, 7 = 0.70, /i p ~ 0.75, 
while /i Q cannot be constrained because of the minimal 
contribution from oblates. Similarly, for the "less lumi- 
nous" sample, fi = 0.45, fi p = 0.70, fi t ,/3 = 0.90 and 
fi t ,-y = 0.70 - 0.75 are derived. 

McMillan et al. (2007) pointed out that the axis ra- 
tios of equal-mass merger remnants can be different for 
differing merging conditions. Then, our result (no signif- 
icant difference between the two samples in terms of the 
axis ratios preferred) could be interpreted as a lack of 
significant difference in the merger history between the 
samples. However, if it takes numerous merging events to 
build elliptical galaxies, the memory of the past merger 
history could easily be buried. The last major/minor 
merger event on the other hand could still be important. 





Fig. 6. — Left: Preferred intrinsic axis ratio within a 1 a range. 
We confirm that the preferred axis ratio of each type converges at 
some axis ratio. Right: Same as top, but for la width of Gaussian 
weight. (T , o-p and at.-y shows good convergence while <Tj a is 
widely accepted. 




Fig. 7. — Apparent ARDs for the two different luminosity sam- 
ples. The solid line and dashed line correspond to "luminous" , and 
"less luminous" samples, respectively. The "luminous" sample is 
rounder than "the less luminous" sample in general. 



6. DEPENDENCE ON ENVIRONMENT 
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Fig. 8. — Same as Fig. [5]but for two different luminosity groups 
(see text). Upper panels show the ARDs for the "luminous" sample 
and lower panels are for the "less luminous" samples. The majority 
of the "luminous" sample are likely triaxial, while a large amount 
of oblate types still exist in the "less luminous" sample. 

Within the context of ACDM scenario, galaxies build 
up hierarchically via galaxy mergers and interactions. In 
this regard, the density-morphology relation (Dressier 
1980) may reveal the significance of environment for 
galaxy formation and evolution. We investigate whether 
the intrinsic shapes of early-types are also related to 
environment. We used Joo Heon Yoon's density pa- 
rameter (p) for the local density of early-type galaxies 
(priv. comm.). Yoon's density parameter (p) is a mea- 
surement of crowdedness counting all the neighboring 
galaxies with a Gaussian-weighting scheme in proportion 
to the distance between the galaxies. Yoon's measure- 
ments are improved over Schawinski et al.'s primarily 
by including more candidate member galaxies assuming 
that early-type galaxies obey optical color-magnitude re- 
lations. Surprisingly, we did not find any notable impact 
of the density parameter on the apparent ARD. Even be- 
tween the two extreme sub-samples (representing fields 
and dense clusters), the apparent ARDs are found to 
share the same parent distribution via a K-S test. It 
should be noted that Ryden (1993) reported a similar 
result, finding no significant difference in the ARD in 
her sample of brightest cluster galaxies (BCGs). 

It is interesting to note that semi-analytic models 
also show consistent results. Khochfar & Silk (2006) 
assert that the stellar properties of merging remnants 
of massive galaxies above the characteristic mass (~ 
3 x 10 10 M Q ; Kauffmann et al. 2003) arc not much dif- 
ferent between the field and cluster environments. Con- 
sidering the brightness of our sample galaxies M r = — 19 



through —23, and assuming a stellar mass-to-light ratio 
5, most of our galaxies would exceed the characteristic 
mass. 

7. LIMITATIONS 

In this section we investigate the limitations of our 
approach by comparing our results with different classi- 
fication schemes. A robust result that does not strongly 
depend on the choice of schemes would ensure the use- 
fulness of the method. Table 1 and Fig. [TU] show the 
three cases defined as Ml, M2, and M3, where M2 is the 
Franx et al. classification scheme we used in this study. 

In Fig. [HI we show the apparent ARDs for the three 
cases assuming a uniform intrinsic ARD. The peak of 
each type can change slightly with the classification 
scheme, but the overall behavior is still retained. In addi- 
tion, the goodness of fit assuming uniformly-distributed 
axis ratios is still poor for other cases {xled ~ 50; see 
§4.1). 

To further constrain the Gaussian parameters, we 
tested the volume-limited sample using three different 
classification schemes. If the dispersion of the preferred 
axis ratio is too large, no fitting distribution can main- 
tain its full information of the assumed distribution. For 
example, a Gaussian distribution of oblate galaxies pa- 
rameterized with [i and a (in the 7 direction) would in- 
evitably have galaxies lying outside the oblate territory 
by the classification scheme (Figure 1): mostly into triax- 
ial. This is particularly so when a is large and for triaxial 
ellipticals due to the type definition. Nevertheless, our 
exercise showed that the preferred axis ratios and type 
weights are robustly derived by assuming Gaussian dis- 
tribution. 

In Table 2, we present the preferred axis ratios (/x) 
corresponding to the most probable case showing reli- 
able convergence and the Gaussian width (cr) for three 
different classification schemes. The results on p and a 
appear to be reasonably consistent between Ml and M2, 
but M3 shows a larger difference. The slight difference 
in the OPT fractions derived is easy to understand. For 
example, the fractions of oblate and prolate in the M3 
case are greater than those in Ml; this is natural because 
there are more possible oblate or prolate configurations 
in M3 than in Ml. In Fig. [lOl we show a mock intrinsic 
ARD for a sample of 4000 galaxies simulated with the 
values in Table 2. We first populate the 4000 x Wtypc 
galaxies for each type using Gaussian distribution. Note 
that oblate and prolate populations have Gaussian distri- 
bution along only one (shortest) axis; oblates (prolates) 
have Gaussianity on 7 (/3). Fig. [TOldemonstrates how the 
intrinsic shapes of galaxies derived change with classifica- 
tion scheme. Once again, Ml and M2 are similar, but not 
M3. This demonstrates the currently unsatisfying situa- 
tion that the results of our OPT analysis can be sensitive 
to the choice of the OPT classification scheme. In this 
sense, it is critical to use a classification scheme that is 
more physically motivated from kinematic requirements. 
If numeric divisions on axis ratios are still the easiest 
scheme, the OPT classification could be determined from 
models. For example, an ensemble of stellar orbits dom- 
inated (by a certain value) by oblate stellar orbits may 
have a characteristic range in c/b/a axis ratios which can 
then serve as the unique criterion for "oblate" . 
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TABLE 1 

Three different classification schemes to test the 
robustness of our opt analysis 



Model 


Oblate 


Prolate 


Triaxial 


Ml 
M2 
M3 


< T < 0.15 
< T < 0.25 
< T < 0.35 


0.85 < T < 1.0 
0.75 < T < 1.0 
0.65 < T < 1.0 


0.15 < T < 0.85 
0.25 < T < 0.75 
0.35 < T < 0.65 


TABLE 2 

THE PREFERRED AXIS RATIOS AND FRACTIONS OF OPT TYPES 
DERIVED 


Model 


/^o & o f-Lp 






Ml 
M2 
M3 


0.44 0.1 0.70 
0.44 0.1 0.72 
0.44 0.1 0.72 


0.1 0.92 (0.3) a 0.74 0.2 
0.1 0.92 (0.25) 6 0.78 0.2 
0.05 0.96 0.1 0.86 0.25 




Model Wo 


W p 


Wt 




Ml 0.24 ± U - 1U 0.15 ±°- uy 
M2 0.29 ±0 - 09 0.26 ±0 11 
M3 0.32 ± 07 0.31 ±0 11 


0.61 ± uli 
0.45 ±0 - 13 
0.37 ± 012 



"poorly constrained 
'poorly constrained 



3.0 C — ' ' ' r- 

- I mT 




0.5 -/ S k- 

~ /• 

0.0 L__, , , i , , , i , , , i , , , . 

0.2 0.4 0.6 0.8 1.0 

Apparent Axis Ratio 

Fig. 9. — Same as Fig. [2]but for different classification schemes. 
The oblate and prolate types are shifted towards the triaxial peak 
as the classification scheme becomes more generous for oblate and 
prolate. 



8. CONCLUSIONS AND DISCUSSION 

We have investigated on the intrinsic shape distri- 
bution of early-type galaxies. We show that the de- 
projection results can be affected by the details of the 
selection criteria of the observed data. We use isopho- 
tal major radius, luminosity, redshift and fracDev com- 
bined in order to minimize sample biases. We have 
constructed a volume-limited sample of 4,994 early-type 
galaxies. In addition, since even the high fracDev cri- 
terion (> 0.95) does not completely remove the contam- 
ination of non-early-type galaxies, we performed visual 
inspection on the galaxies and finally selected 3,922. 

We use the projection probabilities analytically calcu- 
lated to determine the ratios of oblate, prolate and triax- 
ial (OPT) types. We found no combination of randomly- 
distributed OPT types matches the observed ARD! 

We have tested a hypothesis of Gaussian distribution 




b/o(=/S) 

Fig. 10. — A mock intrinsic distribution for a sample of 4000 
galaxies based on the values in Table 2. This figure shows how 
intrinsic ARD changes with classification scheme (see text). 



in axis ratios. The results show excellent fits to the ob- 
served ARD from our volume-limited sample. In partic- 
ular, our results show that triaxial is the most common 
type (W t — 45% ± 12), and they are round in general 
((H,p — 0.92, /j, t}1 = 0.78). On the other hand, oblate 
prefers a more flattened axis ratio (/i = 0.44) and ac- 
counts for approximately 29% of all early- types. The 
prolate show a peak at fi p = 0.72 with a comparable 
fraction to that of oblate (26%). 

Recent numerical N-body simulations suggest that the 
shapes of dark matter halos favor triaxial or prolate over 
oblate (Dubinski & Carlberg 1991, Barnes 1992, Jing k 
Suto 2002, Bailin & Steinmetz 2005, Novak et al. 2006). 
This is consistent with our result at least in the sense that 
triaxial is favored, while a small discrepancy is found in 
the values of the preferred intrinsic axis ratios. Dubinski 
& Carlberg (1991) proposed that dark matter halos are 
triaxial with < j3 >— 0.71 and < 7 >= 0.50, respectively, 
which is somewhat more flattened than we derived. In 
order to reproduce the observed ARD using these pref- 
erences for triaxial type, specifically-designed distribu- 
tions for oblate and prolate types would be needed. This 
is obviously contrived. Interestingly, our results are 
more comparable with the simulation results of Bailin 
& Steinmetz (2005) that suggested 7 = 0.6 ± 0.1 and 
= 0.75 ± 0.15 for dark matter haloes. It should how- 
ever be noted that the mass distribution of halos with 
cooled baryon can be different from those of pure dark 
matter halos (Gnedin et al. 2004). Novak et al. (2006) 
also pointed out the disparity in the halo shape between 
the stellar and dark matter components. If this is true, 
the apparent agreement between the axis ratios derived 
in our study and those from the dark matter simulations 
may not be significant. Further investigations are called 
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for. 

In the context of merging, merger remnants are mainly 
affected not only by initial orientation but also by the 
mass ratio of the merging galaxies (Naab & Burkert 2003; 
Khochfar & Burkert 2006). Bright early-type galaxies 
are often supposed to be formed via numerous merging 
events in the ACDM cosmology. Hence, we should in- 
vestigate all possible parameter space in order to make 
a realistic comparison with observation. In addition, be- 
cause elliptical-elliptical mergers, as well as spiral-spiral 
mergers, can lead to form an elliptical galaxy, the param- 
eter space is large. Furthermore, there appear to be two 
distinct classes of early-type galaxies separated by their 
luminosity (Bender 1988, Kormendy & Bender 1996). 
These studies indicate that luminous galaxies are typ- 
ically supported by anisotropic velocity and have boxy 
isophotcs while faint early-types show disky isophotcs. 
On this basis, we search for a clue on their formation 
history using the OPT parameters we derived. Our re- 
sults indicate that the "luminous" sample (M r < —21.2) 
are mostly triaxial (W t — 67% ± 13) while the "less 
luminous" sample have a large number of oblate types 
(38% ± 8). Interestingly, no clear difference in the pre- 
ferred axis ratios of OPT types between the two samples. 

Khochfar & Burkert (2003; 2005) suggest that the 
origin of the two classes could come from the different 
types of progenitor being merged (c.f., Valluri & Merritt, 
1998). Their semi-analytic models suggest that bright 
early types tend to form from elliptical-elliptical merg- 
ers. Naab et al. (2006) also pointed out that spiral- 
spiral mergers cannot reproduce the observed fraction of 
anisotropic early-types, and that boxy, anisotropic sys- 
tem can form by binary mergers of early-type galaxies. If 
early-type mergers are the main channels to form bright 
anisotropic galaxies, and if our triaxial galaxies corre- 
spond to anisotropic early-types while oblate types are 
rotationally-supported systems, our results show reason- 
able agreement with Naab et al. (2006) in terms of the 
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ratio of anisotropic to the total number of early-type 
galaxies. The fraction of the anisotropic galaxies in the 
"luminous" and the "less luminous" samples are around 
0.8 and 0.4, respectively (Naab et al. 2006), which are 
arguably similar to our results W t — 0.67 ± 0.13 and 
W t = 0.43 ±0.11. 

Since different merging events are proposed to result in 
different configurations for early- type galaxies, it is also 
of interest to investigate whether there is connection be- 
tween environment and galaxy shape. We do not find a 
clear dependence on the local galaxy density. Consider- 
ing that the galaxy number density parameter generally 
represents the dark matter halo potential for the galaxy 
cluster, this may imply that such a cluster-scale environ- 
ment might have little effect to the intrinsic shapes of 
individual galaxies. This is consistent with the results 
from the semi-analytic study of Khochfar & Silk (2006). 

Our analysis has caveats. The quantitative aspects of 
our simulation results depend on the OPT classification 
scheme. Besides, when the fraction of a type (oblate, pro- 
late, or triaxial) is small, our method fails to constrain 
the preferred values of mean axis ratios and Gaussian 
widths. With the improvement on these caveats and de- 
tailed study of N-body simulations, we hope our work can 
provide a simple framework to investigate the formation 
history of early-type galaxies. 
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